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Abstract 



It is demonstrated that the breather solutions recently discovered in the weakly 
coupled topological discrete sine-Gordon system are spectrally unstable. This is in 
contrast with more conventional spatially discrete systems, whose breathers are always 
. spectrally stable at sufficiently weak coupling. 

o 

q ■ 1 Introduction 

There has long been interest in solitonic field theories defined on spacetimes where time 
. is continuous but space is a discrete lattice. The motivation for this scenario varies from 

physical reality (condensed matter and biophysics, see [j^] and references therein), through 
mathematical convenience (integrable quantum field theory ||) to computational necessity 
(when numerically simulating soliton dynamics on a computer). One particularly interesting 
programme, whose motivation is of the last type, is that of Ward. The aim here is to find 
spatially discrete versions of Bogomol'nyi field theories which preserve the Bogomol'nyi 
properties of their continuum counterparts: a topological lower bound on soliton energy 
saturated by solutions of a first order "self-duality" equation, these solutions occuring in 
continuous families, or "moduli spaces." Such systems appear to capture the qualitative 
features of the continuum soliton dynamics even on very coarse lattices, making them ideal 
candidates for computer simulations. The programme is essentially complete for one spatial 
dimension, that is, such a "topological" discretization has been constructed for any kink- 
bearing nonlinear Klein-Gordon model |Q. Results in higher dimensions are more partial 

I IB- 

The most heavily studied example is the topological discrete sine-Gordon (TDSG) sys- 
tem |?J and its variants ||. Here one has an explicit, continuous one-dimensional moduli 
space of static kink solutions, despite the discrete nature of space, and the kink dynamics is 
similar to continuum kink dynamics even when the kink structure covers only 2 or 3 lattice 
sites. This should be contrasted with the conventional discretization, the Frenkel-Kontorova 
model, in which the kink dynamics is extremely complicated || . 

Kinks and antikinks are not the only solitons of the continuum SG equation. It also 
supports breathers, that is, spatially localized, time periodic, oscillatory solutions. It was 
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recently proved |nj that the TDSG system also supports breathers, provided the lattice 
is sufficiently coarse. So the question arises: to what extent are these breathers similar 
to their continuum counterparts? In this note we will demonstrate one crucial difference: 
discrete breathers are spectrally unstable. 



2 Breathers in the TDSG system 

The TDSG system is defined on spacetime Z x E. Space is a regular discrete lattice of 
spacing h € (0, 2] and spatial position will be denoted by n 6 Z. Time is continuous, 
denoted t. The field ip n (t) evolves according to the differential-difference equation 

4 - h? A + h 2 
= Ah 2 cos ^"( sin VVt+i + sinV>n-i) tt^- sin^ n (cos ip n+ i + cos ^ n _i). (1) 

Note that this tends, as h — > 0, to the usual sine-Gordon PDE (in laboratory coordinates 
(x = nh,t)) for (j) = 2ip. Equation (||) supports a continuous translation orbit of static kink 
solutions, 

ijj n = 2tan- 1 e a(n - ft) , b € E (2) 

where a = ln[(2 + h)/(2 — h)]. The conserved energy associated with ([!]), call it E, satisfies 
E > 1 for all ifi satisfying kink boundary conditions {ip n — » 0, n — > — oo; ^ n — > 7r, n — > oo), 
and the static solutions (|2|) saturate this bound. Note that as h — » 2 the kinks tend to step 
functions, which is why one insists h <2. 

Breather solutions to (||) of any period T > 2tt where T/2tt ^ Z, arise as follows. 
Consider the case h = 2. Then (|l]) reduces to 

V^n = -^sinV>n(oOS'0n+l + COS^n-l), (3) 

which supports [0] the "one-site breather" 

*»W = n = 0, ^ 

where is any T-periodic solution of the pendulum equation 

i9 + sin^ = 0. (5) 

By applying an implicit function theorem argument, in the manner of MacKay and Aubry 
[ 11 1 , one obtains the following [1C]: For all T > 2ir, T ^ 2-7rZ, there exists a continuous 



deformation of the one-site breather (Q) through breather solutions of (||) with lattice spacing 
h € [2 — e, 2], where e > is sufficiently small. So breather solutions exist for all h sufficiently 
close to 2, provided the breather period is nonresonant (not a multiple of 2ir). 

By truncating the system to N = 2m + 1 sites (—N<n< N) and applying a Newton- 
Raphson scheme to Pt — Id, where Pt '■ E 2Ar — ► M. 2N is the period T Poincare return map 
for system (|l|), one can construct these breather solutions numerically |0|. One starts at 
h = 2 — Sh, Sh small, with the one-site breather as an initial guess, then interates the NR 
scheme until Pt — Id converges (within some tolerance) to 0. Repeating this for h = 2 — 25h 
with the newly found breather as initial guess, one works piecemeal away from the h = 2 
limit. In this way a portrait of the existence domain of breathers in the (T, h) parameter 
space has been built up [10]. 
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3 Spectral instability 



One says that a period T solution q(t) of a dynamical system on phase space M is spectrally 
stable if the spectrum of the associated Floquet map Fl : T q r \M — ► T q r T \M = T q ^M is 
contained within the closed unit disk DcC. Here Fl = DPt, that is the period T return 
map of the linearized flow about q. Spectral stability is necessary, but not sufficient, 
for linear and hence practical stability. We shall show that the one-site breather (|j) is a 
spectrally unstable solution of (|3|). Since the spectrum of Fl varies continuously under 
constinuous perturbation, it follows that all breathers, for h sufficiently close to 2, are 
unstable. 

Let ip n (t) be a solution of (jl|), and 5ip n (t) be a solution of the linearized flow about 
ip n (t). Then 



+ 



1 Ah 2 
A-h 2 

Ah 2 
A + h 2 

4h 2 



h 2 A + h 2 

SmV>n(smV>n+l + sin^n-l) + 2 COS tp n {cOS V>n+1 + COS^ n _i)]#n 



Ah 2 

COS^ n (cOsVVi+l<5^n+l + COS1p n -l6lp n -i) 

sinV>n(sinV>n+i&/'n+i + sin Ipn-xSlpn-x). 



(6) 



In particular, if h = 2 and ip n (t) is the one-site breather (Q), one finds that 



—5ip n \n\ > 1 

2(1 + cos"d)8ij) n \n\ = 1 
— cos 1? 5ip n n = 0. 



(7) 



Truncating the lattice to size N = 2m + 1 (our result is independent of m), the Floquet 
map takes the initial data of Sip to the final data (at t = T), 



Fl 



' ^-m(0) " 




" <5^_ m (T) ■ 


<ty-m(0) 






#m(0) 


1 — > 




#m(0) 




5lp m (T) 



(8) 



Of course, Fl is linear by linearity of (^). In fact, since (|l|) defines a Hamiltonian flow, Fl 
is symplectic with respect to the natural symplectic structure on M. 2N , that is 



Fl T Q,Fl = Q, 

where is the block diagonal matrix 

S7 = diag( J, . . . , J) , J = 



1 
-1 



0) 



(10) 



Note that Fl is continuously connected through GL(2N, K) to l2jv (by the time evolution 
(||)), so @ implies det Fl = 1. Symplectomorphicity imposes many more constraints on the 
spectrum of Fl [O], but we shall not need these. 
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Since system (□) is decoupled, Fl is block diagonal with one symplectic 2x2 block Fl n 
per lattice site. To demonstrate instability, therefore, it suffices to exhibit one block with 
an eigenvalue outside D = {z £ C : \z\ < 1}. The offending block is Fl\ (or Fl^i): 



Fh 



where y± , 2/2 form a basis of solutions of 



2/1 CO 2/2(0 

yi (O 2/2(0 



(ii) 



y + -(1 + cos^(t))2/ = (12) 

with yi(0) = 2/2(0) = 1, 2/i(0) = 2/2(0) = 0. Note that, as argued above, detF/i = 1, so if 
|tr -PZi I > 2, the eigenvalues are a real conjugate pair {A, 1/A} with |A| > 1, and ^ n (t) is 
spectrally unstable. 

It is numerically trivial to solve ( |i~2| ) (coupled to the pendulum equation (|5|) for i9(t)) over 
[0, T] and hence compute trFZi. The results are shown in figure 1. Clearly trF/i > 2 for 
all T within our range (the limit T — > 00 is numerically inaccessible), so both FZi and Fl-i 
have a real eigenvalue exceeding 1, confirming our claim that the breathers are unstable for 
h close to 2. 

It is possible, of course, that as h decreases substantially below 2, the two offending 
eigenvalues move into D, so that spectral stability is recovered. To eliminate this possibility, 
we have computed the spectrum of Fl at a sample of points within the breather existence 
domain. This domain, whose construction is described in detail in [jlC|l , is depicted in (uj, h) 
space (uj = 2tt/T being breather frequency) in figure 2, with the sampled points marked by 
crosses. Direct computation of the spectrum of Fl is certainly a much cruder approach to 



numerical spectral stability analysis than that of Aubry [13|. However, since we seek only 
to confirm instability, rather than examine bifurcations between stability and instability, a 
simple approach is justified. 

At each point sampled, Fl was constructed by solving (^) approximately using a 4th 
order Runge-Kutta scheme on an 11 site lattice. Note that no decoupling occurs now 
that h 7^ 2. Numerical accuracy of Fl was tested by measuring the deviation of Fl from 
symplectomorphicity by comparing the euclidean norm ||Err|| of 

Err = Fl T QFl — £2, \\M\\ := (J^Afg)^, (13) 

hi 

with I \Fl 1 1 . In each case it was found that the largest eigenvalue of Fl substantially exceeds 
1 . The results are summarized in table ^. 



4 Conclusion 

This paper reinforces the point that breathers in the weakly coupled TDSG system are really 
quite different from continuum sine-Gordon breathers. Their domain of existence does not 
extend to the continuum limit, their initial profiles do not resemble continuum breathers 
sampled on a lattice [Hj, and they are spectrally unstable. In this last property, they are 
also different from breathers constructed in conventional discrete systems (oscillator chains), 
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UJ 


h 


Err 


||FZ|| 


1 A | max 


UJ 


h 


Err 


||M|| 


1 A | max 


0.40 


1.83 


0.0743 


28.6358 


45.2109 


0.65 


1.89 


0.0062 


4.4025 


4.8864 


0.40 


1.90 


0.0854 


22.3127 


30.8483 


0.70 


1.79 


0.0053 


3.9849 


3.9052 


0.45 


1.87 


0.0916 


15.7065 


21.7614 


0.70 


1.92 


0.0004 


3.9330 


3.5739 


0.45 


1.95 


0.0084 


10.6190 


14.7902 


0.75 


1.83 


0.0031 


3.1875 


2.9886 


0.55 


1.85 


0.0011 


7.6868 


10.1095 


0.75 


1.90 


0.0013 


3.0531 


2.8147 


0.60 


1.75 


0.0040 


5.5526 


6.7439 


0.80 


1.88 


0.0008 


2.1615 


2.2982 


0.60 


1.82 


0.0144 


5.5808 


7.3534 


0.80 


1.96 


0.0006 


2.1080 


2.1469 


0.65 


1.78 


0.0040 


4.4930 


4.8864 


0.90 


1.95 


0.0002 


2.6793 


1.4589 



Table 1: The largest eigenvalue of Fl for various breathers. 



which are ubiquitously stable, at least on sufficiently coarse lattices |13J]. We conclude that, 
while the TDSG system captures the dynamics of kinks and antikinks quite faithfully even 
for large h, the same is not true of kink-antikink bound states, that is, breathers. 
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Figure captions 



Figure 1: Plot of trF/i against breather period T. Note that trF/i > 2 for all T. 

Figure 2: Domain of existence of discrete breathers in (u, h) space. The jagged dotted 
curve is the theoretical lower boundary obtained by a phonon resonance argument. The 
solid curve is the actual lower boundary, generated numerically as in |1C]. The small oo 
region (effectively to < 1 /3) is numerically inaccessible. Crosses mark those breathers whose 
spectral instability was confirmed numerically (see table 1). 
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Figure 1 
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Figure 2 
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